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Basic Concepts of the Poisson Process 


11.1.2 Basic Concepts of the Poisson Process 


The Poisson process is one of the most widely-used counting processes. It is usually used in scenarios where we are counting the 
occurrences of certain events that appear to happen at a certain rate, but completely at random (without a certain structure). For 
example, suppose that from historical data, we know that earthquakes occur in a certain area with a rate of 2 per month. Other 
than this information, the timings of earthquakes seem to be completely random. Thus, we conclude that the Poisson process 
might be a good model for earthquakes. In practice, the Poisson process or its extensions have been used to model [24] 

— the number of car accidents at a site or in an area; 

— the location of users in a wireless network; 

— the requests for individual documents on a web server; 

— the outbreak of wars; 

— photons landing on a photodiode. 


Poisson random variable: Here, we briefly review some properties of the Poisson random variable that we have discussed in the 
previous chapters. Remember that a discrete random variable X is said to be a Poisson random variable with parameter fl, shov 
as X ~ Poisson(fl). if its range is Rx = {0,1, 2, 3, . . . }, and its PMF is given by 



for k E Rx 
otherwise 


Px(k) 


Here are some useful facts that we have seen before: 

1. IfX ~ Poisson(p), then EX = /i, and Var(Af) = p. 

2. If Xi ~ Poisson(/J,i ), for i = 1,2, • • ■ , 71, and the Xi's are independent, then 

X\ + X2 + • • • + X n ~ Poisson(/i 1 + H2 + ■ • ■ + p-n)- 

3. The Poisson distribution can be viewed as the limit of binomial distribution. 


Theorem 11.1 

Let Y n ~ Binomial (jl,p = p{n )) . Let fl > 0 be a fixed real number, and lim ra -K» n P = P Then, the PMF of Y n 
converges to a Poisson(p) PMF, as 71 — V 00 . That is, for any k E {0, 1, 2, . . . we have 



Poisson Process as the Limit of a Bernoulli Process: 


Suppose that we would like to model the arrival of events that happen completely at random at a rate A per unit time. Here is one 
way to do this. At time t = 0, we have no arrivals yet, so 1V(0) = 0. We now divide the half-line [0, Oc) to tiny subintervals 
of length d as shown in Figure 11.2. 


0 S 28 3d' 4<5 5d 68 78 88 9(5 10(5 • ■ ■ t 

Figure 11.2 - Dividing the half-line [0, 00 ) to tiny subintervals of length (5. 


t, 


Each subinterval corresponds to a time slot of length 8. Thus, the intervals are (0, (5], (<5, 2d], (2d, 3d], ■ • ■. More generally, the 
kth interval is ((A; — l)d, /cd] . We assume that in each time slot, we toss a coin for which P(H) = p = Ad. If the coin 
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lands heads up, we say that we have an arrival in that subinterval. Otherwise, we say that we have no arrival in that interval. 
Figure 11.3 shows this process. Here, we have an arrival at time t = kS , if the A;th coin flip results in a heads. 
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Figure 11.3 - Poisson process as a limit of a Bernoulli process. 

Now, let N(t ) be defined as the number of arrivals (number of heads) from time 0 to time t. There are n ~ time slots in the 
interval (0, t]. Thus, N(t ) is the number of heads in n coin flips. We conclude that N(t) ~ Binomial(n,p ). Note that 
here p = A<5, so 


np = 


Thus, by Theorem 11.1 . as 5 —V 0, the PMF of N(t) converges to a Poisson distribution with rate A t. More generally, we can 
argue that the number of arrivals in any interval of length T follows a Poisson^Xr) distribution as S —V 0. 

Consider several non-overlapping intervals. The number of arrivals in each interval is determined by the results of the coin flips 
for that interval. Since different coin flips are independent, we conclude that the above counting process has independent 
increments. 

Definition of the Poisson Process: 

The above construction can be made mathematically rigorous. The resulting random process is called a Poisson process with rate 
(or intensity) A. Here is a formal definition of the Poisson process. 

The Poisson Process 


nXS 



A t. 


Let A > 0 be fixed. The counting process {IV(i), t E [0, oo)} is called a Poisson process with rates A if all the following 
conditions hold: 

1. N(0) = 0; 

2. N(t ) has independent increments; 

3. the number of arrivals in any interval of length T > 0 has Poisson(XT ) distribution. 

Note that from the above definition, we conclude that in a Poisson process, the distribution of the number of arrivals in any 
interval depends only on the length of the interval, and not on the exact location of the interval on the real line. Therefore the 
Poisson process has stationary increments. 


Example 11.1 

The number of customers arriving at a grocery store can be modeled by a Poisson process with intensity A = 10 customers per 
hour. 

1. Find the probability that there are 2 customers between 10:00 and 10:20. 

2. Find the probability that there are 3 customers between 10:00 and 10:20 and 7 customers between 10:20 and 11. 
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Solution 


1. Here, A = 10 and the interval between 10:00 and 10:20 has length T = hours. Thus, if X is the number 
of arrivals in that interval, we can write X ~ Poisson{ 10/3). Therefore, 




P(X = 2) = 


(fY 


2 ! 


0.2 


2. Here, we have two non-overlapping intervals I\ =(10:00 a.m., 10:20 a.m.] and I 2 = (10:20 a.m., 11 a.m.]. 
Thus, we can write 


P 3 arrivals in I\ and 7 arrivals in I 2 ) = 


P [ 3 arrivals in I\ ) • P ( 7 arrivals in /; 


Since the lengths of the intervals are T\ = 1/3 and T 2 = 2/3 respectively, we obtain Ari = 10/3 and 
At 2 = 20/3. Thus, we have 

P ^3 arrivals in Iy and 7 arrivals in I 2 

~ 0.0325 


(f) 


(?)’ 


3! 


7! 


Second Definition of the Poisson Process: 


Let N(t ) be a Poisson process with rate A. Consider a very short interval of length A. Then, the number of arrivals in this 
interval has the same distribution as 7V( A) . In particular, we can write 

P{N{ A) = 0) = e“ AA 

A 2 

= 1 — AA H-A 2 — • • • (Taylor Series). 

£ 

Note that if A is small, the terms that include second or higher powers of A are negligible compared to A. We write this as 


P(N(A) = 0) = 1 - AA + o(A) (11.1) 

Here o( A) shows a function that is negligible compared to A, as A —> 0. More precisely, g{ A) = o( A) means that 


lim ^1 = 0. 

A—10 A 

Now, let us look at the probability of having one arrival in an interval of length A. 


P(N( A) = 1 ) = e~ XA XA 

= AA ^1 — AA + — ■ ■ Yj (Taylor Series) 

= AA + A 2 A 2 + ^-A 3 • • 

= AA + o(A). 


We conclude that 
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P(N(A) = 1) = AA + o(A) (11.2) 

Similarly, we can show that 

P(N( A) > 2) = o(A) (11.3) 

In fact, equations 11.1 . 11.2 . and 11.3 give us another way to define a Poisson process. 

The Second Definition of the Poisson Process 


Let A > 0 be fixed. The counting process {N(t), t £ [0, oo)} is called a Poisson process with rate A if all the following 
conditions hold: 

1. N(0) = 0; 

2. N(t ) has independent and stationar y increments 

3. we have 

P(N( A) = 0) = 1 - AA + o( A), 

P(N( A) = 1) = AA + o(A), 

P(N( A) > 2) = o(A). 

We have already shown that any Poisson process satisfies the above definition. To show that the above definition is equivalent to 
our original definition, we also need to show that any process that satisfies the above definition also satisfies the original 
definition. A method to show this is outlined in the End of Chapter Problems. 

Arrival and Interarrival Times: 


Let N(t ) be a Poisson process with rate A. Let X\ be the time of the first arrival. Then, 

P(X i > t) = P(no arrival in (0, t]) 

= e- At . 

We conclude that 

( 1 - e~ xt t > 0 
F Xl (t) = \ 

\ 0 otherwise 

Therefore, X\ ~ Exponential ( A). Let X 2 be the time elapsed between the first and the second arrival (Figure 11.4). 


I X x X 2 X 3 X 4 • • • t 

0 Ti T 2 T 3 T 4 

Figure 11.4 - The random variables X\, X 2 , • • -are called the interarrival times of the counting process N(t ). 

Let S > 0 and t > 0. Note that the two intervals (0, s] and [s, S + t] are independent. We can write 

P{X 2 > t\X\ = s) = P(no arrival in (s, s + t]\Xi = s) 

= P(no arrivals in (s, s + t]) (independent increments) 

= e- At . 

We conclude that X 2 ~ Exponential^ A), and that X 3 and X 2 are independent. The random variables X\, X 2 , ■ ■ ■ are 
called the interarrival times of the counting process N(t). Similarly, we can argue that all X^'s are independent and 
Xi ~ Exponential^ A) fori = 1,2,3, ■ ■ ■. 

Interarrival Times for Poisson Processes 

If N ( t ) is a Poisson process with rate A, then the interarrival times Xi, X 2 , ■ • ■ are independent and 
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Xi ~ Exponential^ A), for i = 1,2,3, • • ■. 

Remember that if X is exponential with parameter A > 0, then X is a memoryless random variable, that is 

P(X > x + a\X > a) = P(X > x), for a, x > 0. 

Thinking of the Poisson process, the memoryless property of the interarrival times is consistent with the independent increment 
property of the Poisson distribution. In some sense, both are implying that the number of arrivals in non-overlapping intervals are 
independent. To better understand this issue, let's look at an example. 


Example 11.2 

Let 7V(f) be a Poisson process with intensity A = 2, and let X\, X 2 , • • ■ be the corresponding interarrival times. 

a. Find the probability that the first arrival occurs after t = 0.5, i.e., P(X\ > 0.5). 
h. Given that we have had no arrivals before t = 1, find P(X\ > 3). 

c. Given that the third arrival occurred at time t = 2, find the probability that the fourth arrival occurs after t = 4. 

d. I start watching the process at time t = 10. Let T be the time of the first arrival that 1 see. In other words, T is the first 
arrival after t = 10. Find ET and Var(T). 

e. I start watching the process at time t = 10. Let T he the time of the first arrival that I see. Find the conditional 
expectation and the conditional variance of T given that I am informed that the last arrival occurred at time t = 9. 

• Solution 

o a. Since X\ ~ Exponent iaZ(2), we can write 

P(X 1 > 0.5) = e -( 2x0 ' 5 ) 

« 0.37 


Another way to solve this is to note that 

P(X 1 > 0.5) = P {no arrivals in (0,0.5]) = e -( 2x0 5 ) « 0.37 

h. We can write 

P(X 1 > 3|Xi > 1) = P(X 1 > 2) (memoryless property) 

= e - 2x2 

« 0.0183 

Another way to solve this is to note that the number of arrivals in (1, 3] is independent of the arrivals before 
t = 1. Thus, 

P(X 1 > 3|Xi > 1) = P(no arrivals in 
= P(no arrivals in 

= e - 2x2 

« 0.0183 

c. The time between the third and the fourth arrival is X 4 ~ 
probability is equal to 

P(X 4 > 2|Xi + X 2 + X 3 = 2) = P(X 4 > 2) (independence of the Xi s) 

= e ~ 2x2 

« 0.0183 


(1.3] | no arrivals in (0,1]) 

(1.3] ) (independent increments) 

Exponential[2). Thus, the desired conditional 


d. When I start watching the process at time t = 10,1 will see a Poisson process. Thus, the time of the first 
arrival from t = 10 is Exponential^ 2). In other words, we can write 


T = 10+ X, 
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where X ~ Exponential^ 2). Thus, 


ET = 


Var(T) = Var(X) 

_ 1 

““ 4' 

e. Arrivals before t = 10 are independent of arrivals after t = 10. Thus, knowing that the last arrival occurred 
at time t = 9 does not impact the distribution of the first arrival after t = 10. Thus, if A is the event that the 
last arrival occurred at t = 9, we can write 


10 + EX 

n 1 21 

10 + - = —, 


E[T\A] = E[T\ 

_ 21 
““ T’ 

Var(T| J 4) = Var(T) 
_ 1 
“ 4' 


Now that we know the distribution of the interarrival times, we can find the distribution of arrival times 

Ti=X u 
T 2 = X 1 +X 2 , 

T 3 = X 1 + X 2 + X 3 , 


More specifically, T n is the sum of n independent Exponential^ A) random variables. In previous chapters we have seen 
that if T n = X\ + X 2 + • • • + X n , where the -X+s are independent Exponential (A) random variables, then 
T n ~ Ganrmria(n , A). This has been shown using MGFs. Note that here n £ N. The Gammaijl , A) is also called 
Erlang distribution, i.e, we can write 


T n ~ Erlangen, A) = Gamma(n , A), for n = 1,2,3, • • ■. 

The PDF ofT„, for n = 1,2, 3, • • •, is given by 

\ nj.n —1 p—\t 

f Tn (t)= - TTT , forf>0. 


(n- !)! 


Remember that if X ~ Exponential(X), then 


m = s , 

Var(X) = +. 


Since T n — X\ + X 2 + ■ • ■ + X n , we conclude that 


E[T n ] = nEX\ = j, 
Var (T n ) = nV ar(X„) = 
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Note that the arrival times are not independent. In particular, we must have T\ < T 2 < T 3 < • • ■. 

Arrival Times for Poisson Processes 


If N(t) is a Poisson process with rate A, then the arrival times T\, T 2 , • • • have Gamma{n , A) distribution. In particular, 
for n = 1,2, 3, • • -, we have 

E[T n \ = and Var(T n ) = 

A A^ 

The above discussion suggests a way to simulate (generate) a Poisson process with rate A. We first generate i.i.d. random 
variables X\, X 2 , X%, • ■ ■, where X t ~ Exponential^ A). Then the arrival times are given by 

T\ = X u 
t 2 = x 1 + x 2 , 
t 3 = x 1 + x 2 + x 3 , 


<— previous 
next —* 
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